NINETEENTH PROGRESS REPORT 


® 7 5 - 


i ft o a i 




on 


CALIBRATION AND EVALUATION OF SKYIAB ALTIMETRY FOR 
GEODETIC DETERMINATION OF THE GEOID (Contract NAS9-13276, 
EPN 440) September 1 to November 30, 1974 

to 


NASA Johnson Space Center 
Principal Investigation Management Office 
Houston, Texas^., 77058 


from 

BATTELLE ' 

Columbus Laboratories 


™ interest of early and wu 0 

semination of Earth Re sources- Survey 
rogram information and without liability 
tor Any use made thereof/' 


December 12, 1974 


Prepared by: S. Gopalapillai and M. Kuhner 

and 

A. G. Mourad (Principal Investigator) 

Z. H. Byrns, Code TF6 - NASA/JSC Technical Monitor 


f7%~) 5- 100 6 4 ) CALIBRATION AND EVALUATION H75-13364 

for SKY.LAB ALTIMETRY FOR GEODETIC 
[DETERMINATION OF THE GEOID Progress 

(Report, 1 Sep. * 30 (Bat telle Columbus Unclas 

(Labs., Ohio . ) 18 p BC $3.25 CSCL 08B G3/43 00064 J 


BATTELLE 

Columbus Laboratories 
505 King Avenue 
Columbus, Ohio 43201 



NINETEENTH PROGRESS REPORT 


on 

CALIBRATION AND EVALUATION OF SKYLAB ALTIMETRY FOR 
GEODETIC DETERMINATION OF THE GEOID (Contract NAS9-13276, 
EPN 440) September 1 to November 30, 1974 

to 

NASA Johnson Space Center 
Principal Inves tigation Management Office 
Houston, Texas 77058 

from 

BATTELLE 

Columbus Laboratories 

Pi^epg.x'ed. by: S. Gopalapillai and. M. Kuhner 

and 

A, G. Mourad (Principal Investigator) 

Z. H. Byrns, Code TF6 - NASA/JSC Technical Monitor 
December 12, 1974 


PR0GR.ES S 

During this period the tapes containing the modified SKYBET data 
have been received. The main effort during this period was aimed at 

(1) matching the altimetry data from ECT S071-1 tapes (received 
in August) with satellite ephemeris data from SKYBET tapes with respect to 
the time of observation, 

(2) modifying the computer program used in the simulation studies 
to make it more general, efficient and economical with respect to computer 
time to suit the data on tapes and to enable the solution of a larger system 
of equations simultaneously, and 

(3) testing this program with the real data. 
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The details of matching the altimetry data with the ephemeris 
data are presented in Appendix A. The reduced algorithm for modifying 
the existing computer programs are presented in Appendix B. The data and 
documents received, during this period, are reported in Appendix C. 

DATA PROCESSING RESULTS 


In the earlier analysis, only one bias term was considered for 
all the data from various submodes. The obvious discrepencies in the 
heights between different submodes had been manually accounted for. The 
modified algorithm assumes a different bias term for each submode so that 
these terms can absorb the height discrepencies in the submodes. This 
assumption is well supported by the results obtained from the real data 
analysis. These results include: 

(1) the recovered bias terms for each submode are distinctly 
different. The magnitude of the differences agrees with the height 
discrepencies in the observed data, 

(2) the insensitivity of the bias terms to the a priori weights 

of the observations indicate the stability and efficiency of the mathematical 
model used for evaluating them, and 

(3) the accuracy of estimates obtained from the analysis are 
highly correlated with the number of available ground truth data points. 

PROBLEMS AND RECOMMENDATIONS 


There are no significant problems or recommendations to be made 


for this period. 
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NEXT PERIOD 

The major efforts during the next reporting period include: 

(1) Continuation of the processing and analysis of data 
from the remainder of EREP Pass No. 9, 

(2) Evaluation of the geoidal profile obtained for Pass 
No. 9 with the ground truth profile, and 

(3) Initiation of the processing and analysis of data 
from Passes No. 4, 6, and 7. 

TRAVEL 

No plans for travel are anticipated at this time. 



APPENDIX A 


MATCHING OF ALTIMETRY DATA WITH SKYBET DATA 


For each altimeter range to be corrected, the precise location 
in space of the satellite is required. Therefore, altimetry data points 
must be combined with ephemeris orbit data points from the SKYBET tapes. 
Unfortunately the times for which these data are recorded are different 
for the two types of tapes. On the SKYBET tapes the time interval between 
data points is exactly 1/8 second. On the altimetry tapes the time interval 
is 1.04/8 seconds. The result is that the two sets of data wander 
slowly in and out of phase. The data needed from the SKYBET tapes 
are the satellite latitude, longitude, height above reference ellipsoid 
and the earth-fixed geocentric X, Y, and 2 coordinates. The first three 
of these parameters do not change significantly in 1/8 second. However, 
the X, Y, and Z coordinates which must be known to better than 1 meter do 
change by several meters in 1/8 second. Therefore, a simple interpolation 
scheme was used to compute the X, Y, and Z coordinates to better than one 
meter for the times at which the altimeter ranges are given. 

For each time, t, on the SKYBET tape not only X, Y, and Z 

v • * 

coordinates are given but also the velocity components X, Y, and Z. 
Therefore, at a later time t+At the approximation, 

X( t+A t) - X(t) + At X(t) 

Y(t+At) = Y(t) + At Y(t) 

Z(t-hH) = Z(t) + At Z(t) 

can be used. It can be shown that if no forces other than gravity affect 
the satellite, X, Y, and Z are very nearly constant over 1/8 second so that 
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this approximation is valid (the errors should be of the order of 0.06 
meters or so). To test the interpolation scheme, the X, Y, and Z 
coordinates were computed for At = 1/8 second and the results were compared 
with those given on the SKYBET tapes. The maximum error was 0.095 meter, 
which is of the same order of magnitude as estimated from theory. Therefore, 
this interpolation scheme was incorporated into a computer code which 
combines the necessary SKYBET data with the necessary altimetry data and 
writes them onto a single tape for later processing. 



APPENDIX B 


MODIFIED ALGORITHM FOR COMPUTER PROCESSING 
OF ALTIMETRY DATA 

1.0 CONDITION EQUATION OF INTRINSIC PARAMETERS 

cl 

The adjusted value R., corresponding to the measured altimeter 

range, R°, is intrinsically related to (a) the geocentric coordinates, 

v . y . 7 . of the satellite at the instant of measurement, (b) the geoidal 

sr» si’ si’ 

undulation, , at the satellite subpoint referred to a given reference 

r 

ellipsoid, and (c) the biases in all the measurement systems involved. 

There are two types of condition equations considered in this 
investigation, depending on how the system biases are modeled. In the case 
where the system biases are considered proportional to the measured ranges, 
the condition equation is given by 

R a ( 1-hAf ) + N a - D. =0 (1) 

i xi 

where 


D i “ F (x si’ Y si’ z si’ a ’ e) 


(la) 


being the height of the satellite above the reference ellipsoid given as 
a function of the geocentric coordinates of the satellite at the instant 
of observation , and the parameters (a, e) defining the size and shape, 
respectively, of the assumed reference ellipsoid. R^Af is system bias 

for this observation. Hie second type of equation is 


R a + Af + N a - D. = 0 
i li 


( 2 ) 


where the bias Af is assumed to be constant independent of the range. 
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The first form was used in the. studies reported in the previous 
progress reports. However, the vai'iation in R^ is very small compared to 
R. itself so that R.Af remains constant for all practical purposes. 

Preliminary results of the numerical comparisons of these two models support 
this contention. Therefore, the model given in equation (2) will be used 
throughout the rest of this investigation. 

2.0 THE OBSERVATION EQUATION 

In order to determine the bias, Af, the other quantities in 
equation (2) are needed: is the measured altimetry range; is the 

geoid undulation available as part of the ground truth data, and is 
assumed to be known for each of the measured ranges. Both R_^ and Ik are 
associated with random errors which have to be determined in a least 
squares adjustment. 

The observation equation for the model described in equation 
(2) is given by 

V. + Af + AN. + W. - 0 (3) 

i ii 

where V and AN are residuals of R. and N. respectively and 

ii ii 

W. = R° + N° - D. (4) 

1111 

where the superscript o refer to the observed values of the quantities 

R and N . Equation (3) refers to only one altimetry range. There will be 
i i 

as many equations of this type as there are altimetry observations for 
which the ground truth data is available. The system of all these equations 
can be written in matrix notation as follows. 
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V + AX + W = 0 (5) 

where 
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Preliminary examinations of the altimetry range ovservations 

for different submodes indicate different biases for the observations in 

th 

each submode. In this case equation (3) for the i observation which 

til 

corresponds to the k submode will be of the form 


V. + Af + AN. + W. - 0 
1 k 1 1 


( 10 ) 


which will result in the A matrix of the form 


Af , Af rt Af, An AN, 


A = 


0 0 


AN 


n 


( 11 ) 


while the structure of other matrices will remain unchanged, 


3.0 MATRIX SOLUTION 


If the weight matrices for the observations R and N are P and 
respectively, their least squares solution in matrix form is given by 


- (A 1 PA + P .) _1 A 1 PW 

v X 

(12) 

- AX-W 

(13) 


where the superscript 1 referes to the transpose of the matrix. 


The variance covariance matrix, Z , for the vector X is given by 


Z 

x 


a 2 (A 1 pa + P ) 
o x 


with being 


the variance of unit weight given by 


( 14 ) 
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0 (X 1 P X + V X PV) 

a 2 = 5 

o df 

where, df is the degree of freedom. 

4.0 SPECIAL CHARACTERISTICS OF THE MATRICES INVOLVED 


4.1 Weight Matrices 


( 15 ) 


The observations R and the a priori estimates for N are assumed 

to be independent. The biases corresponding to each of the submodes are 

also assumed to be independent of each other. Under these assumptions, the 

matrices P and P are diagonal. If the P matrix is partitioned along a 
x x 

line separating the biases from the undulations then, P can be written 

X 

in the form 

(16) 

Further, if the accuracy of all the observations in each of the 
groups R, N, and Af is assumed to be equal then, 




A 


l 

I 


I 


(18) 
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where is the submatrix of A containing the columns corresponding to 
the bias terms Af . Similarly X is partitioned into components X and 

K 

X^ such that 


(19) 


The special structures of the matrices, as described above, are 
taken advantage of in the numerical evaluation of equation (12) - (15). 


5.0 REDUCED ALGORITHM FOR DIGITAL COMPUTER EVALUATION 


The nature of the weight matrices assumed and the structure of 
the design matrix on either side of the partition make it possible to 
simplify the equations given in equations (12) - (15) for computer coding 
so that this program could handle the data more efficiently and economically. 

The partition of matrices enables the solution of the normal 
equation (12) to be sequential, i.e., to solve for x^, and then for x^. 
Equation (12) can be rewritten as 

X - - N _1 U (20) 

wi th 

N = (A 1 PA + P v .) 

( 21 ) 

and U = A PW 


Using the partition approach, the submatrices of N and U will be 


— 



— [ 


N 11 

! N 12 


pA^i + px r L f 

pAj 

N 21 

1 N 22 
1 JJ 

1 

p\ ! 

(p+px 2 ) . I 


N 


( 22 ) 



B-7 


U 


Now, let 


r " 


. 1 „ 

u i 


PA x W 

— 



U 2 


pW 

_ _ 



Q 11 ! Q 12 


N n ' 

N 12 


= 

4 - — i 

Q 21 ] Q 22 


l> . 

N 22 1 


-1 


Then (Faddeev and Faddeeva, 1963) 

Q 


11 


(N u - N 1Z N-l 


Q 12 = " Q 11 N 12 N 22 

Q 22 = (I - nj 2 Q 12 ) 


and the solutions X^ and X 2 are given by 
X, = -Q. 


11 


u i ’ N i2 n 22 u ; 


= - N 22 N U X 1 • K 22 U 2 


(23) 


(24) 

(25) 

(26) 

(27) 

(28) 


With the subraatrices of N defined as in equation (22) in terms of the 
weight matrices and the design matrix A, it can be shown that 


11 


P • P x 2 a 1 A + ox 
• A i i p i.i 


(29) 


p+px. 


It should be noted that since A^ is diagonal, Q;q is also diagonal, 


W22 U 2 


P* P x 2 

— A7W 

p+px^ 1 


and U 


(30) 
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Then, from equations (27) and (28) 


X 1 ” - Q u u 


p+px. 


[ A i x i* 


( 31 ) 

(32) 


The weight coefficients matrix for the geoid undulations N (Ground Truth) 
is given by which is given by 


•22 


p+px. 


X + 


p+px. 


Vu A I 


(33) 


One of the advantageous features in this system is that both the matrices 
(Q and N 97 ) whose inverses are required are diagonal. This enables one 
to solve a system (of normal equations) of any size with relatively small 
computer storage requirement. 

Once the bias for each submode is determined, the geoid undulation 
for any other point along the pass is determined from equation (2) as 

N = D - R - Af (34) 


6.0 TEST RESULTS AND DISCUSSIONS 

Using the algorithm presented in the last section, the data 
corresponding to the first 3 submodes in mode 5 of EREP pass 9. (Start 
Time 163:13:1:37.181, Stop Time 163:13:4:26.701) was evaluated. 

The input data are: 

(a) The altimetry observations at 3 seconds interval. 
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(b) The SKYBET ephemeris data corresponding to these 
altimetry data. 

(c) Geoid undulations from Marsh-Vincent geoid map. 

(d) The weight matrices corresponding to the altimetry 

observation (p), bias (px^) , and to the geoid 
undulations • 

The 5 seconds interval for altifnetry observations is selected 
simply because this is the least interval for which the satellite sub-points 
can be plotted on the available Marsh-Vincent geoid map. 

The reference ellipsoid used is the one used by Marsh and 
Vincent (1974) where it is defined by 

semi-major axis a = 6378142 m 

flattening f = 1/298.255 

Five test runs were made varying the relative weights p» px^, 

2 

and P x 2 * The resu l ts are presented in Table 1, where a is a posteriori 
variance of unit weight. 


TABLE 1 



Test 1 

Test 2 

Test 3 

Test 4 

Test 5 

p 

1/1.0 

1/9.0 

1/9.0 

1/4.0 

1/1.0 

px l 

1/100,0 

1/100.0 

1/03 

1/co 

1/co 

px 2 

1/0.01 

1/4.0 

1/4.0 

1/4.0 

1/4.0 


-20.75 

-20.15 

o 

oO 

0 

CM 

1 

-20.80 

-20.80 

A £ 2 

-7.26 

-7.22 

-7.27 

-7.27 

-7.27 

Af 3 

2.81 

2.77 

2.81 

2.81 

2.81 

a 

9.15 

0.92 

0.85 

1.05 

1.35 


o 
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The a priori estimate assumed for the variance of unit weight 


is unity. Comparing this value to the a posteriori variance, it is felt 
that the relative weights between the altimetry observations and the 
ground truth are equal (tests 2, 3, and 4). However, the estimates 
obtained for the bias terms are significantly insensitive to the relative 


weights used. 

The a posteriori standard deviations obtained for the bias terms 


are 


Bias 


Af. 


Af, 


Af , 


Standard deviation 
(in meters) 

1.50 

0.66 

0.95 


No. of points of Ground Truth 
4 

21 

10 


These values indicate that the standard deviation is inversely 
proportional to the square root of the number of ground truth points 
available for the observations in each submode. 


APPENDIX C 


REPORTS AND DATA RECEIVED 


Title 


. Identification 
Number 


(1) SKYLAB PROGRAM EARTH RESOURCES EXPERIMENT PACKAGE 

July 10, 1974 Contract NAS8-24000 

Amendmen t J SC - 14 S 

Sensor Performance Report, Vol IV (S193 R/S) 
(Engineering Baseline, SL2 and SL3 evaluation) 

(2) NASA TECHNICAL MEMORANDUM 

June, 1974 NASA TM X-58122 

(3) SKYLAB EREP 

S191 INFRARED SPECTROMETER 
DATA ACQUISITION CAMERA SCENE LIST FOR 
SL-2 , SL-3 , and SL-4 


(4) 


(5) 


( 6 ) 


(7) 


SKYLAB 4 S190A 461636 461536 PI 

4X Transparencies - 1 each pos. 


Mag : 

52 

FRAMES: 083/090 

Mag : 

70 

FRAMES: 194/207 . 

Mag : 

A4 

FRAMES: 366/462 

EREP 


Tape No. 

4 


V14612 

6 


V03546 

7 


V05607 

9 


V03248 

54 


V03284 

85 


V03208 

97 


V09770 

SKYLAB 4 

October 

21 

2 x Prints - 1 each 
Mag 90 

057/60 61 thru 62 blank 

SKYLAB PROGRAM 

EARTH RESOURCES EXPERIMENT PACKAGE 

Sep temb 

er 6 , 

1974 Contract NAS8-24000 


Amendment JSC-14S 

Sensor Performance Report, Volume 1 (S190A) 
(Engineering Baseline, SL2, SL3 and SL4 Evaluation) 


MSC-055 29 
(SL3) 


JSC-08461 

440 


W06079 


MSC-05528 


\ 


\ 

\ 


No . of 
Copies 

1 

1 

1 

1 

7 

4 

1 


C-2 


Title 


Identification 

Number 


No. of 


(8) SKYIAB PROGRAM EARTH RESOURCES EXPERIMENT PACKAGE MSC-05528 

September 6, 1974 Contract NAS8-24000 

Amendment JSC-14S 

Sensor Performance Report, Volume II (S191) 

(Engineering Baseline, SL2, SL3 and SL4 Evaluation) 

(9) REQUIREMENTS FOR EREP ELECTRONICS SENSORS PHO-TR524 

October 18, 1974 Rev. A, Ch. 2 

(10) SKYLAB PROGRAM EARTH RESOURCES EXPERIMENT PACKAGE MSC-05528 

September 6, 1974 Contract NAS8-24000 

Amendment JSC-14S 

Sensor Performance Report, Volume V (S193 ALT) 

(Engineering Baseline, SL2, SL3 and SL4 Evaluation) 

(11) SKYLAB PROGRAM EREP MSC-05528 

September 6, 1974 Contract NAS8-24000 

Amendment JSC-14S 

Sensor Performance Report, Volume V (S192) 

(Engineering Baseline, SL2, SL3 and SL4 Evaluation) 

(12) EARTH RESOURCES EXPERIMENT PACKAGE (EREP) MSC-07744 

EXPERIMENT CALIBRATION DATA 


I 

3 

1 

1 

I 



